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introduction 


- This  memorandum  describes  a ray  tracing  program  which  has  been  written 
for  the  Datatron  205  to  permit  the  computation  of  corrected  spreading  loss 
for  convergence  zone  rays  in  a layered  ocean.  This  program,  called  Ray 
Tracing  Program  VI,  differs  from  Program  II  in  two  major  respects:  first, 
it  includes  a computation  of  the  diffraction  loss  correction  factor, 
needed  in  the  vicinity  of  each  caustic  to  give  a finite  value  of  spreading 
loss;  and  second,  it  limits  printouts  to  those  ranges  at  which  the  ray 
passes  through  the  receiver  depth.  Rays  are  delineated  by  their  values 
of  vertex  velocity,  V*,  rather  than  vertical  angles  at.  the  source,  and  four 
types  of  ray  paths  may  be  computed,  with  a variable  number  of  range  cycles 
between  source  and  receiver. 


THEORY 


The  subject  of  corrections  to  be  used  in  the  vicinity  of  caustics  was 
introduced  by  Marsh  in  Study  H of  reference  (a).  It  was  shown  that  the 
complete  expression  for  spreading  loss  in  this  case  should  contain  a mul- 
tiplicative term,  |3Y«C)|&,  called  the  first-order  diffraction  correction. 
This  rather  complicated  correction  term  is  a function  of  frequency  and 
modifies  the  classical  spreading  loss  equation,  so  that  the  intensity  of 
a ray  at  a caustic  remains  finite,  even  though  dR/dVx  approaches  zero. 


The  variable  in  the  correction  factor  is  given  by 
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where  cC  ~ correction  factor  variable 
R = range 

SQ  = vertical  angle  of  the  ray  at  the  source 
/e  = sound  velocity  at  the  source 
f = frequency  in  cps. 

Since  Program  VI  doas  not  include  a choice  of  frequency  at  which  the 
correction  is  to  be  made,  the  output  contains  the  correction  factor  , 

which  can  then  be  converted  to  cC  . Figure  1 gives  the  corresponding 
value  of  diffraction  loss  in  db,  which,  when  added  to  the  classically 
derived  spreading  loss  given  in  the  printout,  yields  the  corrected  value 
of  spreading  loss  for  the  ray  considered. 

RAY  TRACING  PROGRAM  VI 

The  Datatron  20b  program  to  be  described  traces  rays  from  a source 
depth  to  a receiver  depth,  separatee  by  a prescribed  number  of  deep  re- 
fractions and  shallow  refractions  or  reflections,  and  prints  information 
on  range,  travel  time,  spreading  loss,  and  the  diffraction  loss  correction 
factor,  . 


INPUT  REQUIREMENTS 


1.  A layered  velocity-depth  profile  (0n>  Vn)»  with  the  depth  in  fed  and 
the  velocity  in  ft/sec.  The  program  will  handle  a maximum  of  forty  layers. 


2.  The  depth  of  the  source  (D6),  which  must  corresnonc  to  a depth  on  the 
velocity-depth  profile  of  paragraph  1. 

3.  The  depth  of  the  receiver  (Dr),  which  must  correspond  to  a depth  on 
the  velocity-depth  orofile  of  paragraph  1. 


4,  The  vertex  velocity  of  the  ray  (Vx)»  which  ittermines  the  vertical 
angle  at  any  point  along  the  path.  Several  vertex  velocities  can  be  com- 
puted in  the  same  problem  by  specifying  the  first  and  last  rays  and  the 
increment  between  adjacent  rays. 


f>.  The  ray  is  terminated  after  a specified  number  of  "range  cycles’  (n). 
The  first  oassage  between  Ds  and  Dr  depends  upon  the  "ray  type"  chosen 
(see  paragraph  6),  and  then  n complete  traverses  through  maximum  ano  mini- 
mum depths  are  added,  ending  at  Dr. 
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b.  There  are  four  ray  types,  of  which  any  combination  can  be  chosen  for 
computation,  as  follows* 

Type  i travels  down  from  Ds  and  up  to  Dr 
Type  2 travels  up  from  Ds  and  up  to  Dr 
Type  3 travels  down  from  D6  and  down  to  Dr 

Type  4 travels  up  from  Ds  and  down  to  Dr 

7.  Travel  time  (T)  is  an  optional  feature  in  the  program,  which  can  be 
included  in  the  computations  by  suitable  coding  of  the  input  data. 

OUTPUT 

The  output  of  Program  VI  may  be  a printed  list  or  a set  of  punched 
cards  The  format  of  tie  printed  output  is  as  follows! 

Line  1 - 

Word  1 - Depth  of  source,  Dt. 

Word  2 - Depth  of  receiver,  Dr. 

Line  2 - 

Word  1 - Vertex  velocity,  Vx- 

Word  2 - Depth  at  which  ray  vertexes  (deepest  point). 

Lins  3 - 

Word  1 - Ray  type  (l,  2,  3,  cr  -1). 

Line  4 - 

Word  1 - 0 cycles  (at  this  point,  the  ray  is  traced  from  Ds  to  Dr 
via  the  path  specified  by  the  ray  type). 

Word  2 - Horizontal  range,  R,  in  yards. 

Word  3 - dR/dVx,  Ir,  yds/ft/sec.  sec  reference  (b)). 
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Word  4 - Vx  • d^R/dVj^,  in  yds/ft/sec » 

Word  5 - Spreading  loss  in  db. 

Word  6 - , the  diffraction  loss  correction  factor. 

Word  7 - Travel  time,  I,  in  seconds  {if  desired). 

Lines  5 through  n - 

Word  1 - The  number  of  complete  range  cycles  added  to  the  above. 

Words  2 through  7 - The  corresponding  values  of  R,  dR/dVx,  Vx’d^R/dV;^, 

NSpr,  <vVk  , and  T which  result  from  the  combination 
of  the  basic  ray  path  (Line  4)  followed  by  the 
specified  number  of  complete  range  cycles. 

Lines  3 through  n will  be  repeated  for  each  ray  type  requested. 

The  format  listed  above  is  repeated  from  Line  2 on  for  each  additional 
value  of  Vx  requested. 

The  following  exceptions  to  the  above  Printout  are  made; 

1.  If  a requested  vertex  velocity  is  less  than  or  equal  to  the  velocity  at 
the  source,  Vs,  then  Vx  and  Vs  will  be  printed  on  Line  2 and  the  program 
will  advance  to  the  next  value  of  Vx  requested. 

2.  for  any  ray  for  which  Dr  does  not  fall  within  the  upper  and  lower  depth 
limits  of  the  ray,  there  will  be  no  special  printout.  After  Line  2 is 
printed,  the  program  will  advance  to  the  next  value  of  Vx  requested. 

3.  All  rays  must  vertex  below  the  source,  for  any  ray  which  bounces  off 
the  bottom,  the  Vx  anc  a word  of  all  93s  will  be  printed  on  Line  2. 

EXAMPLE 

To  test  the  ability  of  Program  VI  to  produce  acceptable  results,  an 
example  was  selected  for  which  previous  computations  involving  diffraction 
corrections  had  been  made.  The  velocity-depth  profile  given  in  Table  I of 
reference  (c),  as  shewn  in  Figure  2,  was  used  as  the  input  to  Program  VI, 
and  the  results  for  all  four  ray  types  in  the  first  convergence  zone  were 
obtained  for  a source  depth  of  50  feet  and  a receiver  depth  of  100  feet. 

The  uncorrected  spreading  loss  fox-  each  ray  type  is  shown  in  Figure  3.  The 
diffraction  corrections  for  a frequency  of  530  cps  were  applied  to  these 


rays,  resulting  in  the  corrected  values  of  spreading  loss  shown  in  Figure  4. 
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This  figure  agrees  quite  well  with  the  results  given  in  Figure  2 of 
reference  ( c ) . To  show  the  effect  of  combining  ail  of  the  available  rays, 
assumption  of  random  phase  between  arrivals  was  made,  and  Figure  5 re- 
sulted. The  solid  curve  was  derived  from  the  data  of  Figure  3,  and  the 
dashed  curve  from  the  data  of  Figure  4.  Figure  5,  again,  agrees  satis- 
factorily with  Figure  3 of  reference  (c),  thus  confirming  the  validity  of 
the  program  and  the  correction  method  described. 

CONCLUSIONS 

Ray  Tracing  Program  VI  can  be  used  in  deriving  the  corrected  spreading 
loss  for  rays  in  the  vicinity  of  caustics.  In  its  present  form,  this  pro- 
gram suffers  from  the  necessity  to  choose  values  of  Vx  in  the  input  direc- 
tions, so  that  it  is  often  necessary  to  make  two  or  three  runs  with  the 
same  velocity  profile  before  adequate  coverage  of  a convergence  zone  results. 
This  feature  would  be  eliminated  by  a computed  choice  of  sequential  /alues 
of  Vx,  which  could  be  based  on  the  desired  maximum  spacings  of  points  in 
the  output  plot.  An  attempt  will  be  made  in  the  near  future  to  derive 
additional  program  steps  to  give  this  most  desirable  feature, 
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